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We study the performance of a Wolff-type embedding algorithm for RP^ cr-models. We find that the algorithm in 
which we update the embedded Ising model a la Swendsen-Wang has critical slowing-down as 2^^. fii 1. If instead we 
update the Ising spins with a perfect algorithm which at every iteration produces a new independent configuration, 
we obtain ~ 0. This shows that the Ising embedding encodes well the collective modes of the system, and that the 
behaviour of the first algorithm is connected to the poor performance of the Swendsen-Wang algorithm in dealing 
with a frustrated Ising model. 



In recent years there has been a lot of work 
in devising new algorithms which, by taking into 
proper account the collective modes of the the- 
ory, are able to eliminate or at least to reduce 
critical slowing-down. 

For 0{N) cr-models, an extremely efficient al- 
gorithm was proposed three years ago by Wolff 
. In two dimensions, numerical tests of the dy- 
namic critical behaviour show the complete or 
almost complete absence of critical slowing-down 
(i.e. z < 0.1) 

The extraordinary efficiency of this algorithm 
has spurred many attempts to find generaliza- 
tions to cr-models taking values in manifolds 
other than spheres. However, last year we 
presented a heuristic argument whose conclusion 
was: a necessary condition for a WolfF-type em- 
bedding algorithm to work well (even with per- 



fect updating of the induced Ising spins) is that 
the manifold is a sphere, a real projective space, 
or a discrete quotient of products of such spaces 

i- 

Let us briefly review the general principles of 
Wolff-type embedding algorithms |^,|| . Consider 
a general cr-model taking values in a Riemanian 
manifold M , with Hamiltonian of the form 

H{{a}) ^ (3j2^(''-'''yy (1) 

Then the algorithm is defined by a collection of 
energy-preserving maps T, and gives rise to the 
induced Ising Hamiltonian 

- ^ hxyiex - ey) + const (2) 

i^y) 
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where {e} are Ising spins and 

Jxy = ^[E(Ta.^,(Ty) + E{a^,Tay) - 2E(G^,ay)\ 

^^y = ^[E{Ta^,ay) - E{a^,Tay)] (3) 

In practice an iteration of the algorithm works 
as foUows: 

(i) Choose a map T in the given family according 
to a given distribution. 

(ii) Initialize all Ising spins e^; = 1. 

(iii) Update the embedded Ising model. 

(iv) Set (Jx = Tax where ex ~ —1. 

In step (iii) one can use any valid algorithm for 
simulating the Ising model (2). We will consider 
two different choices: 

a) The practical algorithm where step (iii) 
consists of one standard (full-lattice) Swendsen- 
Wang update. 

b) The idealized algorithm where at step (iii) 
we generate a new configuration of Ising spins, 
independent of the old one. This is achieved in 
practice by performing at every iteration Nhu 
Swendsen-Wang updates, where N^it is chosen 
so large that the autocorrelation times of the var- 
ious observables are independent of Nhn within 
error bars. 

The idealized algorithm allows us to under- 
stand how well the embedding succeeds in em- 
bodying the important large-scale collective modes 
of the cr-model. A bad performance of the ideal- 
ized algorithm means that in the u-model there 
are other important excitations which are not 
captured by the embedding. By contrast, a poor 
performance of the practical algorithm might be 
due solely to the bad performance of the algo- 
rithm used in updating the Ising spins. 

What we have defined is a generalization of 
Wolff's algorithm for 0{N) models, and we claim 
that it can work well only in a few cases. The 
reason for this is that in order to perform well 



the algorithm must do a good job in handling the 
collective modes of the theory, which certainly 
include long- wavelength spin- waves. In order to 
treat these modes well, we argue that the set of 
links for which Jij « must disconnect the x- 
space into two or more regions. It follows that 
the embedding map must have the codimension- 
1 property: the fixed-point manifold of the map T 
must have codimension 1. Differential geometry 
can then be used to prove that the only manifolds 
which satisfy this requirement are or RP^ 
(and discrete quotients of products thereof). 

Let us notice that our heuristic argument 
gives a necessary condition for the idealized (and 
hence also the practical) algorithm to beat criti- 
cal slowing-down, but it does not guarantee that 
either the idealized or the practical algorithm 
will in fact perform well. For this reason we 
have decided to study the two-dimensional RP^ 
model. 

The real projective space RP^~^ is by def- 
inition the sphere S^^^ with antipodal points 
identified, i.e. RP^-^ = S'^-^/Z^. The most 
convenient approach is to consider spins taking 
values on the sphere S*^"^, subject to the con- 
dition that the Hamiltonian and all physical ob- 
servables must be invariant under the Z2 local 
gauge transformations Ux — > rfxCfx with rjx = ±1- 
The simplest lattice Hamiltonian for this model 
is therefore 

^(W) = -f E(^--^-+A.)' (4) 

The continuum limit of this model is not at all 
clear. In the formal continuum limit a — > 0, 
the Hamiltonian becomes that of the contin- 
uum 0{N) non- linear cr-model. In order to ex- 
plain why in the continuum limit the theory does 
not have the Z2 gauge invariance, it has been 
suggested |^ that at a finite value of the cou- 
pling the system undergoes a phase transition 



S. Caracciolo et al. / Wolff-type algorithms for RP^ a-models 



3 



which gives rise to a condensation of the vortices. 
However, the presence of this phase transition is 
rather controversial (see |^,0,||j9) and references 
therein). We do not have yet much to add to 
this point, and in the following we will address 
the problem of the dynamical behaviour of the 
algorithm. 

The algorithm is defined by the same embed- 
ding used by Wolff for 0{N) cr-models: the in- 
duced Hamiltonian is given by (2) with h^y — 
and 

Jxy Picr^ ■ <7y)i<Ja; ■ r){ay ■ r) , (5) 

where — cr^ — ((T-r)r. Let us notice that, when 
> 3, the induced Hamiltonian is frustrated. 
Let us first discuss the behaviour of the prac- 
tical algorithm. We have measured the energy, 
the tensor susceptibility xt and the correlation 
length in the tensor channel ^ for both RP^ and 
RP^ on lattices of dimension L = 32,64, 128 (a 
detailed discussion of the simulation is given in 
0). 

A finite-size scaling analysis of L~^^Tint,x '^^i'" 
sus ^/L shows that the points are well fitted us- 
ing 



0.9 ± 0.3 for i?p2 



1.1 ±0.3 for i?p3 
while a similar analysis for the energy gives 



Zint,E — 



0.2 ± 0.3 for 



(6) 



(7) 



0.2 ± 0.3 for 
This means that the practical algorithm, though 
providing a significant improvement over local al- 
gorithms, still suffers from strong critical slowing- 
down. At this point, however, it is not clear what 
is the cause of this behavior: are there other exci- 
tations in the model which are not well encoded 
in the embedding, or is the critical slowing-down 
due instead to the Swendsen-Wang subroutine 
which is unable to simulate efficiently a frus- 
trated Ising model? 



To answer this question we have studied the 
idealized algorithm for RP^ on lattices with L = 
32,64. We have found that the dynamic critical 
exponent for the susceptibility is now 



0.1 ±0.3 



(8) 



Critical slowing-down is thus nearly eliminated! 
We conclude that the embedding encodes well 
the collective modes of the RP^ model, and that 
the failure of the practical version must be as- 
cribed to the Ising subroutine. 
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